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Outage-Efficient  Strategies  for  Multiuser 
MIMO  Networks  with  Channel  Distribution 
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Department  of  Electrical  and  Computer  Engineering 
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Abstract —  In  this  work,  we  examine  single  user  and  mul¬ 
tiuser  Multiple-Input  Multiple- Output  (MIMO)  beamform¬ 
ing  networks  with  Channel  Distribution  Information  (CDI). 
Since  CDI  changes  infrequently  compared  to  Channel  State 
Information  (CSI),  algorithms  based  on  CDI  can  achieve 
significant  savings  in  feedback  compared  to  algorithms 
based  on  CSI.  With  CDI,  we  can  only  guarantee  quality  of 
service  for  a  specified  outage  probability  in  the  network. 
Assuming  correlated  Rayleigh  fading  on  all  the  links,  we 
derive  a  closed-form  expression  for  the  outage  probability. 
Then,  using  this  expression,  we  derive  algorithms  for 
joint  transmit/receive  beamforming  and  power  control  to 
minimize  the  weighted  sum  power  in  the  network  while 
guaranteeing  these  outage  probabilities.  For  both  single- 
user  and  multiuser  MIMO  scenarios,  we  present  optimal 
algorithms  under  the  Kronecker  model  assumption,  and  we 
present  near-optimal  algorithms  assuming  general  correla¬ 
tion  structures  on  the  links.  We  then  show  that  using  these 
algorithms  based  on  CDI,  if  we  are  willing  to  accept  given 
outages  on  the  links,  we  can  achieve  comparable  power 
usage  in  the  network  relative  to  algorithms  based  on  CSI. 

Index  Terms — covariance  feedback,  beamformers,  mul¬ 
tiuser  MIMO,  outage  probability,  Channel  State  Infor¬ 
mation,  Channel  Distribution  Information;  EDICS:  MSP- 
CAPC,  MSP-MULT,  WIN-ADHC,  WIN-PHYL 

I.  Introduction 

Multiuser  multiple-input  multiple-output  (MU- 
MIMO)  networks  have  generated  much  research  interest 
in  recent  years.  As  the  community  has  come  to  have 
a  greater  understanding  of  MIMO  point-to-point  links, 
recent  work  has  involved  looking  at  how  multiple 
antennas  can  be  utilized  to  reduce  interference  in 
multiuser  scenarios.  Some  useful  surveys  and  texts  on 
the  topic  are  [l]-[4].  To  achieve  maximum  throughput 
and  reliability  in  the  MU-MIMO  network,  all  the  nodes 

This  work  was  sponsored  by  the  U.  S.  Army  Research  Office  under 
the  Multi-University  Research  Initiative  (MURI)  Grant  #  W911NF-04- 
1-0224  and  the  USARAC  Grant  #  W91  INF-09- 1-0494. 


need  to  have  perfect  Channel  State  Information  (CSI) 
of  all  the  links  in  the  network.  With  this  knowledge, 
utilizing  intelligent  beamforming,  interference  in  the 
system  can  be  minimized  and  optimal  rates  can  be 
achieved. 

This  assumption  may  be  feasible  in  small  networks 
where  the  channels  change  slowly.  The  nodes  then  have 
time  to  estimate  their  channel  and  feed  back  their  CSI 
to  the  rest  of  the  nodes  before  transmitting  their  data. 
However,  in  reality,  many  networks  are  large,  mobile, 
and  operate  in  dynamic  environments.  Thus,  the  channels 
change  too  fast  for  the  CSI  to  be  fed  back  to  all  the 
nodes-by  the  time  all  the  CSI  has  reached  every  node, 
the  information  is  already  outdated.  Thus,  much  work 
has  focused  on  finding  transmission  schemes  that  do 
not  rely  on  full  CSI.  The  simplest  transmission  scheme 
assumes  no  knowledge  about  the  channel  structure  and 
employs  equal  power  allocation  on  all  the  antennas  [1]. 
This  approach  has  the  advantage  of  having  no  feedback 
overhead,  but  yields  lower  throughput  and  gives  no 
insight  into  the  transmit  power  required  for  reliable 
multiuser  communication  as  compared  to  schemes  that 
utilize  feedback.  Another  transmission  scheme  uses  lim¬ 
ited  CSI,  where  nodes  feed  back  their  CSI  only  to 
nearby  nodes  [7].  This  approach  assumes  that  the  nearby 
nodes  dominate  the  interference,  and  should  therefore 
yield  near-optimal  performance.  However,  even  in  point- 
to-point  systems,  very  low  channel  coherence  times 
make  full  CSI  feedback  infeasible.  Many  schemes  thus 
utilize  quantized  forms  of  CSI  to  reduce  feedback  [5]  [6]. 
Still,  these  approaches  suffer  from  having  to  feed  back 
information  every  time  the  channel  changes. 

Another  approach,  then,  is  to  utilize  the  channel 
statistics,  or  Channel  Distribution  Information  (CDI), 
to  enhance  communication.  Since  it  takes  into  account 
the  randomness  in  the  channel,  CDI  is  more  robust  to 
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small  channel  coherence  times  and  is  thus  valid  for  much 
longer  than  CSI.  In  addition,  given  the  trend  towards 
more  location-based  services,  statistical  data  based  on 
node  location  can  also  be  collected  and  stored  a  priori. 
Such  location-based  data  eliminates  the  need  for  real¬ 
time  channel  feedback.  CDI  has  been  combined  with 
various  CSI  schemes  to  improve  performance  and  reduce 
feedback  in  [8]-[10],  In  [8],  CDI  is  used  to  aid  in  SINR 
feedback,  and  in  [10],  CDI  is  used  in  combination  with 
channel  norm  feedback  to  exploit  multiuser  diversity. 
In  [9],  CDI  is  utilized  for  more  efficient  beamformer 
codebook  design.  While  these  schemes  use  CDI  to 
reduce  CSI  feedback,  channel  information  must  still  be 
fed  back  every  time  the  channel  changes.  Thus,  much 
previous  work  has  focused  on  utilizing  CDI  without  the 
aid  of  CSI. 

Previous  work  using  only  CDI  focuses  on  maximizing 
ergodic  capacity.  In  [11],  optimal  transmission  strategies 
are  given  for  single-user  multiple-input  single -output 
(SU-MISO)  channels  when  either  the  mean  channel 
information  or  covariance  channel  information  is  known. 
The  capacity-optimal  input  covariance  for  the  single-user 
MIMO  (SU-MIMO)  channel  with  mean  and  covariance 
channel  information  is  presented  in  [12].  This  result  is 
extended  to  the  MIMO  multiple  access  channel  in  [13]. 
Suboptimal  solutions  for  the  MIMO  broadcast  channel 
are  given  in  [14].  These  works  focus  on  average  capacity, 
but  this  metric  relies  on  averaging  over  good  and  bad 
channels  to  achieve  the  desired  throughput.  Thus,  packet 
delay  is  ignored  in  such  schemes,  and  average  capacity 
analysis  is  not  applicable  to  delay-sensitive  applications. 

This  work  seeks  to  account  for  this  delay  problem 
by  looking  at  the  outage  of  the  links  in  the  network. 
Due  to  the  randomness  in  the  channel,  reliable  trans¬ 
mission  cannot  be  achieved  all  the  time  using  just 
CDI.  Consequently,  in  this  work  we  guarantee  a  cer¬ 
tain  signal-to-interference-plus-noise  ratio  (SINR)  with 
a  specified  probability  on  all  the  links.  Specifically, 
this  work  looks  to  jointly  optimize  the  power  alloca¬ 
tion  and  transmit/receive  beamformers  to  minimize  the 
power  used  in  the  network  while  meeting  these  outage 
SINR  constraints.  While  many  works  in  the  literature 
assume  the  channels  experience  independent  identically- 
distributed  (i.i.d.)  Rayleigh  fading  across  the  transmit 
and  receive  antennas  for  analytic  convenience,  much 
work  has  been  done  showing  the  insufficiency  of  this 
model  [15]-[18],  Therefore,  this  work  assumes  all  links 
in  the  network  experience  correlated  Rayleigh  fading. 
Under  this  assumption,  the  expression  for  the  outage 
probability  on  each  link  is  derived.  This  expression  is 
then  applied  to  the  SU-MIMO  channel,  and  an  algorithm 
for  joint  transmit  and  receive  beamforming  and  power 
control  is  presented.  We  then  study  the  MU-MIMO 


system,  where  algorithms  for  joint  transmit  and  receive 
beamforming  and  power  control  are  discussed.  Under  the 
Kronecker  model  assumption  [15],  an  optimal  algorithm 
is  given,  and  for  general  correlation  structures  a  close- 
to-optimal  algorithm  is  derived.  Then,  using  the  derived 
CDI  algorithms,  we  show  that  while  we  accept  some 
loss  on  the  links  in  the  network,  we  can  achieve  power 
consumption  levels  attained  by  CSI  schemes. 

The  paper  is  organized  as  follows:  in  Section  II,  the 
problem  formulation  is  given  and  discussed,  and  the 
expression  for  outage  probability  is  derived.  Section  III 
discusses  the  SU-MIMO  problem.  Section  III- A  finds 
the  closed-form  expression  for  the  optimal  power  control 
for  a  fixed  set  of  beamformers.  Section  III-B  derives  the 
optimal  beamformers  for  the  Kronecker  model,  and  III- 
C  derives  an  algorithm  for  beamforming  under  general 
Rayleigh  fading.  Section  IV  discusses  the  MU-MIMO 
problem.  First,  IV-A  derives  an  optimal  power  control 
algorithm  with  a  fixed  set  of  beamformers.  Section  IV-B 
then  derives  the  optimal  beamformers  for  the  Kronecker 
model,  and  Section  IV-C  discusses  a  good  suboptimal 
solution  for  the  beamformers  assuming  general  Rayleigh 
fading.  We  compare  our  derived  CDI  algorithms  to  CSI 
algorithms  in  Section  V,  and  we  summarize  our  results 
in  Section  VI. 

In  this  work,  the  following  notation  is  used:  italicized 
letters  indicate  scalars  (e.g.  pi ),  lower-case  bold  letters 
indicate  vectors  (e.g.  v),  and  upper-case  bold  letters 
indicate  matrices  (e.g.  A).  Furthermore,  (»)H  indicates 
the  Hermitian  operator,  and  (•)*  indicates  the  conjugate 
operator. 

II.  Problem  Formulation 
A.  System  Model 

This  work  considers  time-varying  MIMO  channels  for 
many  users  in  a  network.  Consider  a  MIMO  network 
with  L  transmit-receive  pairs.  At  link  /,  the  transmitter 
sends  the  symbol  s;(f)  to  the  receiver.  The  transmitter 
uses  unit-norm  beamforming  vector  v/  (t)  to  precode  the 
signal,  and  transmits  with  power  piit).  The  receiver  em¬ 
ploys  the  linear  unit-norm  beamformer  u/ (t)  to  combine 
the  signal.  The  channel  from  transmitter  i  to  receiver  I 
is  given  by  H u(t).  The  noise  A)(i)  is  distributed  as  a 
complex  circular  Gaussian,  and  represents  the  combined 
noise  after  applying  the  receive  beamforming  vector  to 
the  incoming  signal.  The  Ith  received  signal  is  thus  given 
by 

n{t)  =  y/pi(t)[u^  (t)Hu(t)vi(t)\si(t)  +  Ni(t) 

L 

+ v/piw[uf  (f)Hti(f)vi(f)]si(f) 
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In  schemes  that  use  perfect  CSI,  a  block-fading  model 
is  assumed,  so  the  channel  stays  constant  over  each 
block.  Then,  for  notational  convenience,  the  time  vari¬ 
able  will  be  dropped  for  the  channel,  power  allocations, 
and  beamformers.  To  further  simplify  notation,  define 
Gu  =  |ufHh;v!;|2  as  the  beamforming  channel  gain 
from  the  transmitter  on  link  i  to  the  receiver  at  link  l 
and  rr2V(  as  the  noise  power  for  the  Ith  link.  Then,  under 
this  model,  the  SINR  1’/  on  each  link  can  be  shown  to 
be 


PiGu 

PiGli  +  0JV, 


(1) 


If  perfect  CSI  is  available,  to  ensure  a  reliable  link 
is  available  to  all  nodes  in  the  network,  each  link  has 
an  SINR  constraint:  F/  must  be  greater  than  a  threshold 
7 1.  The  goal  is  then  to  minimize  the  power  consumed  by 
the  network  while  meeting  all  the  SINR  constraints.  The 
cost  function  considered  in  this  work  is  the  weighted  sum 
power.  In  this  setup,  each  link  l  in  the  network  incurs 
some  cost  wi  >  0  to  transmit  across  its  link.  An  example 
of  a  network  with  varying  power  costs  on  the  links  are 
networks  with  varying  battery  life  at  the  transmitters. 
For  minimizing  non-weighted  sum  power,  w/  =  1  for 
l  =  l,...,L. 

To  compact  notation,  define  the  weighting  and  power 
vectors  as  w  =  {wi, . . .  ,Wl}  and  p  =  {pi, . . .  ,Pl}, 
respectively.  The  beamforming  matrices  are  defined  as 
U  =  {u1;...,uL}  and  V  =  {v1;...,vL}.  The  op- 
timization  problem  for  having  perfect  CSI,  formulated 
and  discussed  in  [7],  can  then  be  stated  as  follows: 


T 

mm  w  p 

p>o,u,v 

s.t.  =  (2) 

For  a  fixed  set  of  channel  matrices  H/,;’s,  this  problem 
can  be  solved  and  will  give  a  set  of  power  allocations 
and  beamformers  for  the  transmitters  and  receivers  in 
the  network.  Then,  for  every  change  in  the  H;,'s,  all  the 
transmit  and  receive  beamformers  must  be  updated,  and 
the  power  allocation  scheme  changes.  In  many  networks, 
the  feedback  required  for  these  changes  in  the  channel 
is  unrealistic  due  to  rapidly-varying  CSI.  Thus,  this 
work  considers  a  network  where  instantenous  CSI  is 
unavailable,  but  CDI  is  available.  When  only  CDI  is 
known,  the  exact  Hj,’s  are  not  known-instead,  they  are 
assumed  to  be  a  random  variable  drawn  from  a  complex- 
normal  distribution: 


sider  the  case  where  the  channel  varies,  but  the  statistics 
stay  constant.  Under  this  model,  the  expression  for  SINR 
given  in  (1)  becomes  a  random  variable  since  it  depends 
on  the  channel.  The  constraints  in  (2)  can  then  no  longer 
be  written  as  the  SINR  on  link  l  always  exceeding  some 
threshold  7;-since  the  SINR  is  now  random,  it  will 
drop  below  7 1  with  some  probability.  Therefore,  these 
absolute  constraints  change  to  outage  constraints,  and 
links  are  allowed  to  have  an  SINR  below  their  thresholds 
for  specified  probabilities.  Mathematically,  the  constraint 
on  link  l  in  (2)  becomes  Pr(T;  <  7 ;)  <  a;,  where  a; 
is  the  probability  that  the  link  is  in  outage.  The  main 
optimization  problem  using  CDI  can  then  be  formulated: 

T 

mm  w  p 

p>0,U,V 

s.t.  Pr(rI<7l)<a«,Z  =  l,...,L  (3) 

In  this  formulation,  the  constraints  ensure  that  the  Ith 
link  will  be  active  with  probability  (1  —  a;).  In  order  to 
solve  this  optimization  problem,  the  first  step  is  to  derive 
a  closed-form  expression  for  the  outage  probability, 
Pr(Ti  <  7 1). 


B.  Derivation  of  Outage  Probability 

Two  steps  are  used  in  deriving  this  closed-form  ex¬ 
pression  of  the  outage  probability.  The  first  is  to  express 
U  as  a  ratio  of  an  exponential  random  variable  to  a 
weighted  sum  of  exponential  random  variables  plus  some 
constant.  The  second  step  is  to  express  Pr(F;  <  77)  in 
terms  of  the  moment-generating  function  for  a  sum  of 
weighted  exponential  random  variables  to  obtain  the  final 
result. 

Without  loss  of  generality,  the  expression  for  outage 
probability  of  the  first  user  will  be  shown  (treating 
all  other  users  as  interference).  Furthermore,  for  con¬ 
venience,  the  user  subscript  l  will  be  dropped  (Gu 
becomes  G 1,  0%  becomes  and  so  on).  Thus,  the 
expression  for  SINR  being  considered  is 


PiGr 

PiGi  +  a2N 


(4) 


The  following  Lemma  will  also  be  used  in  the  proof  of 
the  theorem: 

Lemma  1:  If  u  and  v  are  unit-norm  vectors  and 
H  is  complex  circular  Gaussian  matrix  distributed  as 
vec(H)  ~  CN( 0,S),  then 


vec(Hw)  ~  CN(0,  Ek) 

The  channel  covariance  matrices,  given  by  the 
comprise  the  CDI  of  the  network.  This  work  will  con- 


uhHv| 


(v*  (g>  u)ffE(v*  (gi  u) 


xi 


(5) 


The  proof  is  given  in  Appendix  I.  The  theorem  can  then 
be  stated  as  follows: 
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Theorem  1:  In  a  MU-MIMO  network  where  all  links 
experience  correlated  Rayleigh  fading  and  the  transmit¬ 
ter  and  receiver  both  employ  linear  beamforming,  the 
expression  for  outage  probability  for  the  SINR  for  user 
1  is  given  by 

Pout  =  Pr(r  <  7)  =  + 

L\  v  c1Plj 

(6) 

where 

Ci  =  (v*  ®  Ui)HSj(v*  (8)  ui) 

Proof:  Step  1\  Get  pout  into  the  following  form: 

Pout  =  Pr(r  <  7)  =  Pr(  ^  <  7),  (7) 

r  +  cjz 


pOO 

Pout=  (1  -  e~i(v+a2))fY{y)dy 
Jo 

nOO 

=  1  -  e_  20'2  /  e~ivfY{y)dy  (10) 

Jo 

Here,  the  key  observation  is  that  the  expression  in  the 
integral  in  Eqn.  (10)  is  the  moment-generating  function 
of  Y,  denoted  ipY(t),  evaluated  at  —3).  Since  Y  is  a  sum 
of  independent  random  variables,  the  moment-generating 
function  of  Y  is  the  product  of  the  moment-generating 
functions  of  each  element  in  the  sum.  The  moment¬ 
generating  function  of  an  exponential  random  variable 
is  given  as 

ipZi{t)  =  E[etZi]  =  (1  —  2f)_1 

The  moment-generating  function  of  Y  is  then  given  as 


where  X  is  an  exponential  random  variable  and  Y 
is  a  weighted  sum  of  independent  exponential  random 
variables.  Applying  Lemma  1  to  the  G;  terms  in  the 
expression  for  T  in  Eqn.  (4)  gives 


r  =  cipiZi 

Yh= 2  ciPiZi  +  aN  ’ 
Zi  ~  xl ,Z'iS  i.i.d. 


(8) 


Then,  dividing  the  top  and  bottom  of  the  right  hand  side 
of  Eqn.  (8)  by  cipi  gives 


r  = 


X 

Y  +  o2’ 


=  E[etY]  =  E[e(^kiZi]  =  E\etk'zf 

i—2 

L  L 

=  Yl.i’Ziikt)  =  J|(i  -  2  hty1 

i—2  i—2 

Then,  evaluating  i/)Y(t)  at  t  =  —  substituting  the  result 
into  Eqn.  (10),  and  substituting  gives  the 

result: 


Pout  =  1  -  e  2<7  Y1 
i= 2 


where  X  =  Zlt  Y  =  2  kiZi,  h  =  and  cr2  = 

.  With  these  substitutions,  pout  has  the  desired  form 
of  Eqn.  (7).  This  concludes  Step  1. 

Step  2:  This  step  expresses  pout  in  terms  of  the 
moment-generating  function  of  Y  to  obtain  the  final 
result.  The  analysis  here  uses  techniques  similar  to  work 
done  in  [19], [20],  Consider 


Pout  =  Pr(vf:  2  <  7)  =  Prpf  <  7 {Y  +  cr2)) 

Pr(2f  <  7 (y  +  a2))fY(y)  dy  (9) 

Here,  observe  that  Pr(X  <  7 (y  +  cr2))  is  the  cumulative 
distribution  function  (CDF)  of  the  exponential  distribu¬ 
tion  (A  =  7)  evaluated  at  7 (y  +  o2),  which  gives 

Pr(X  <  7 (y  +  a2))  =  1  —  e_2^+<T  1 

Substituting  this  expression  into  Eqn.  (9)  yields 


C.  Observations  on  the  Optimization  Problem 

Using  the  expression  for  outage  probability  above  and 

substituting  cr;  =  the  optimization  problem  in  (3) 
can  be  written  as 


mm  w  p 

p>0,U,V 


_  2 L  •/vz 
s.t.  1  —  e  2  a“Pi 


i^l 


n  1+ 


■a- 


CliPi 
'  cupi 


-1 


^  Oil  )  l  \  ,  .  .  .  y  L 

(ID 


To  get  the  constraints  in  (11)  into  a  more  convenient 
form,  define 


5;(7;,p,U;,  V) 

A 


'll  "*1 

=  e  2  cnPl 


n 


i^l 

=  (l-«z)-1 


CliPi  \ 
CllPl  ) 
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Then,  with  some  manipulation  on  the  constraints,  (11) 
can  be  rewritten  as 

T 

mm  w  p 

p>0,U,V 

s.t.  gii'ji,  p,  u;,  V)  <  /3i,l  =  1, . . . ,  L  (12) 

The  form  presented  in  (12)  will  serve  as  the  main 
optimization  problem  used  for  analysis  in  the  rest  of 
this  work.  To  gain  more  insight  into  this  problem,  two 
lemmas  are  given  that  will  aid  in  the  algorithms  and 
manipulations  for  future  sections.  These  lemmas  are 
extensions  of  the  work  in  [7]  to  this  problem.  The  proofs 
are  given  in  Appendix  I. 

Lemma  2:  Holding  the  beamformers  U  and  V  con¬ 
stant,  the  optimization  problem  in  (12)  has  the  following 
3  properties: 

1)  If  it  exists,  the  optimum  solution  p*  is  unique. 

2)  All  the  outage  constraints  are  active  (i.e.  hold  with 
equality)  for  p  =  p*. 

3)  p*  is  unaffected  by  the  choice  of  w. 

Intuitively,  all  the  constraints  are  active  since  if  any 

of  them  held  with  strict  inequality,  the  objective  func¬ 
tion  could  be  reduced  by  making  the  constraint  hold 
with  equality.  Also,  since  all  the  constraints  hold  with 
equality,  there  are  L  equations  (the  constraints)  and 
L  unknowns  (the  pi’s).  Thus,  the  constraints  uniquely 
determine  the  optimal  solution  p*,  and  so  the  weighting 
vector  has  no  effect  on  the  solution.  It  follows  from 
Lemma  2  that  for  fixed  transmit  and  receive  beam- 
formers,  solving  the  constraints  with  equality  will  yield 
the  optimal  power  allocation  scheme,  and  the  objective 
function  in  (12)  can  be  ignored.  The  next  related  Lemma, 
analogous  to  [7],  will  aid  in  gaining  insight  into  subop- 
timal  solutions  in  Section  IV-C. 

Lemma  3:  Consider  the  following  optimization  prob¬ 
lem: 


T 

max  w  p 
p>0,U,V 

s.t.  5/(7z,p,uz,  V)  >  pi,l  =  1, . . .  ,L  (13) 

This  problem  has  the  same  solution  as  the  optimization 
problem  in  (12). 

Intuitively,  the  optimization  problem  in  (13)  is  the 
same  as  the  problem  in  (12)  since  once  again,  all  the 
constraints  will  be  active,  so  the  constraints  uniquely 
determine  p*. 

This  foundation  aids  in  the  development  of  algorithms 
to  find  solutions  for  the  optimization  problem  in  (12). 
First,  the  simpler  single  user  MIMO  problem  is  consid¬ 
ered. 


III.  The  Single  User  MIMO  Problem 
A.  Optimal  Power  Control 

Despite  the  rich  literature  in  single  user  MIMO,  the 
problem  where  both  the  transmitter  and  receiver  only 
know  statistical  information  about  the  channel  and  want 
to  achieve  a  certain  SINR  threshold  with  some  prob¬ 
ability  has  not,  to  our  knowledge,  been  studied.  Here, 
there  is  no  interference  in  the  system,  and  the  weighted 
cost  function  does  not  apply  since  the  problem  requires 
minimizing  over  one  variable.  Thus,  the  optimization 
problem  for  this  special  case  is  as  follows: 

min  P 

P>  0,u,v 

s.t  exp  (  t  s)  <  P  (14) 

Although  there  are  no  interference  terms,  it  can  be  seen 
that  the  proof  for  Lemma  2  holds,  so  the  constraint  in 
(14)  is  satisfied  by  equality.  Solving  the  constraint  for  P 
gives 


2(v*  ®  u)ffS(v*  ®  u)log/3 

The  main  challenge  in  this  problem  is  then  to  find  the 
optimal  beamformers  u  and  v.  From  Eqn.  (15),  it  can  be 
seen  that  P  will  be  minimized  when  (v*®u)ff  S(v*®u) 
is  maximized.  Thus,  the  following  optimization  problem 
follows: 

max(v*  ®  u)'ffS(v*  ®  u) 

11, V 

S.t.  1 1  uj  |  2  =  1,  1 1 V 1 1 2  =  1  (16) 

Solving  for  the  beamformers  in  (16)  is  equivalent  to 
solving  for  the  beamformers  in  (14).  First,  the  Kronecker 
model  will  be  considered,  where  a  closed-form  solution 
can  be  obtained.  Then,  the  general  correlated  Rayleigh 
fading  case  is  considered,  where  an  iterative  algorithm 
to  find  the  beamformers  is  derived. 

B.  Optimal  Beamforming  for  the  Kronecker  Model 

The  Kronecker  model  assumes  that  the  spatial  covari¬ 
ance  matrix  can  be  broken  into  the  spatial  correlation  at 
each  link  end,  so  the  covariance  matrix  is  given  as 

£  =  ££®£ji  (17) 

This  model  is  widely  used,  and  its  validity  is  discussed 
and  verified  in  [15]  [23]  [24] .  To  obtain  the  beamforming 
optimization  problem  for  the  Kronecker  model,  substi¬ 
tute  the  expression  for  £  given  in  Eqn.  (17)  into  (16). 
The  objective  function  in  (16)  then  becomes 
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(v*  <8>  u)ff£(v*  8  u)  =  (v*  8  u)ff(Sy  8  S^)(v*  8  u) 
=(vT£?v*)(uff£fiu)  =  (vff£Tv)(uH£Ku) 

The  second  equality  follows  from  the  property  of  Kro- 
necker  products  that  states 

(A  8  B)(C  8  D)  =  (AC  8  BD),  (18) 

and  the  last  equality  follows  from  aT  =  a,  where  a  is  a 
scalar  (so  vtE£v*  =  (vtE£v*)t  =  vff£Tv).  With 
these  manipulations,  the  maximization  problem  in  (16) 
then  becomes 

max(vffSj’v)(uffI]flu) 

U,V 

s.t.  ||u||2  =  1,  ||v||2  =  1  (19) 

Note  that  u  and  v  are  now  separated,  so  two  separate 
maximization  problems  arise  whose  solutions  are  well 
known:  u  is  the  (normalized)  eigenvector  that  corre¬ 
sponds  to  the  largest  eigenvalue  of  £#,  and  v  is  the 
(normalized)  eigenvector  that  corresponds  to  the  largest 
eigenvalue  of  E t-  It  is  worth  noting  that  if  S t  and  E  R 
are  set  to  the  identity  matrices  of  appropriate  dimension, 
then  E  =  I  and  the  problem  simplifies  to  the  special  case 
of  i.i.d.  Rayleigh  fading.  In  this  case,  any  set  of  transmit 
and  receive  beamformers  are  optimal,  since  all  vectors 
are  eigenvectors  of  I,  and  the  eigenvalues  are  all  1. 

After  they  are  calculated,  the  solutions  for  u  and  v  can 
be  substituted  back  into  the  constraint  of  (14)  to  solve 
for  P.  Define  a!/P  and  as  the  largest  eigenvalues  of 
E t  and  E r,  respectively.  Then  Popt  is  given  by 


(v*  8  u)H E(v*  8  u)  =  vt£uv*  =  uH £vu 

It  can  be  shown  that  £u  and  £v  are  positive  semidef- 
inite  matrices  since  for  a  positive  semidefinite  matrix 
E,  AHEA  is  also  positive  semidefinite  [25].  Now,  £u 
can  be  used  to  solve  for  the  optimal  transmit  beam- 
former  for  a  fixed  receive  beamformer,  and  £v  can  be 
used  to  solve  for  the  optimal  receive  beamformer  for 
a  fixed  transmit  beamformer.  In  this  way,  an  iterative 
algorithm  can  be  derived  to  obtain  both  the  transmit 
and  receive  beamformers.  First,  if  v  is  fixed  such  that 
Ev  8  0,  with  at  least  1  nonzero  eigenvalue  (random 
initialization  should  suffice),  the  optimal  solution  for 
u  will  be  the  eigenvector  corresponding  to  the  largest 
eigenvalue  of  £v.  Then,  if  u  is  then  fixed,  the  optimal 
solution  for  v  will  be  the  complex  conjugate  of  the 
eigenvector  corresponding  to  the  largest  eigenvalue  of 
£u.  This  algorithm  can  alternately  solve  for  u  while 
fixing  v  and  then  solve  for  v  while  fixing  u.  The 
optimization  problem  in  (16)  is  not  jointly  convex  in 
u  and  v,  so  this  algorithm  will  converge  monotonically 
to  a  locally  optimal  solution.  To  initialize  the  algorithm, 
the  Kronecker  model  approximation  can  be  used  on  the 
covariance  matrix,  and  the  optimal  beamformers  from 
this  approximation  can  be  used.  While  this  solution  is 
not  proven  to  be  globally  optimal,  it  is  shown  in  Section 
V  to  perform  well. 

Once  the  algorithm  to  obtain  the  transmit  and  receive 
beamformers  is  applied,  a  solution  for  P  can  be  obtained 
using  Eqn.  (15): 


p  _  Wn 
opt  ~  ~\(D\(i) 


2At  A#  log/? 


(20) 


Popt 


2«pt  ®  uopt)H'Z(Kpt  ®  U opt)  log/? 


(23) 


C.  Beamforming  for  General  Rayleigh  Fading 

Now  consider  the  general  Rayleigh  fading  case.  Using 
the  property  of  Kronecker  products  given  in  Eqn.  (18), 
observe  that 

v*  8  u  =  (v*  8  Ijvh)(1  8  u)  =  (v*  8  Iath)u  (21) 
v*  8  u  =  (Ijvt  8  u)(v*  8  1)  =  (Ijvt  ®  u)v*  (22) 

Then,  define 

Sv  =  (v*  8  IjvK)ffE(v*  8  Ijvh) 

Su  =  (1nt  ®  u)hT,(1Nt  8  u) 

From  Eqns.  (21)  and  (22),  it  follows  that 


IV.  The  Multiuser  MIMO  Problem 

A.  Optimal  Power  Control 

Now  return  to  the  multiuser  problem  in  Eqn.  (12). 
This  problem  is  more  difficult  to  solve  than  the  single- 
user  case  due  to  interference  from  other  users.  In  order 
to  make  the  problem  tractable,  first  examine  the  problem 
when  both  the  transmit  and  receive  beamformers  are 
fixed.  With  this  assumption,  a  power  control  algorithm 
solving  the  constraint  equations  in  (12)  with  equality  is 
derived.  Algorithms  for  the  beamformers  are  then  studied 
in  the  following  sections. 

To  derive  an  optimal  power  control  algorithm,  start  by 
taking  the  logarithm  of  the  constraint  equations  in  (12): 
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—  Nl  t — r  / 

10g(5j(7J>P)U«.  v))  =  log  e  2  cnPi  I  1  +  7; 


=  7^at, 

ZcuPl 

<  log  Pi 


EM1 


Multiplying  both  sides  by 


log  A 


i^l 


results  in 


7i<4, 


Pi 


2 cu  log  Pi  log  Pi 

Now  define 


Elog  1 


i^l 


■7 1- 


CliPi 
'  cupi 


<Pi  (24) 


^(P) 

I(P) 


Wn, 

2  cu  log  Pi 


Pi 

log  Pi 


EloS 

i^l 


The  function  I(p)  is  a  standard  interference  function. 
For  the  definition  of  standard  interference  functions  and 
the  proof,  see  Appendix  II.  A  key  property  of  standard 
interference  functions  is  that  they  satisfy  I(p)  <  p. 
Thus,  using  this  function  and  a  starting  power  vector 
p,  the  update  equation  for  the  algorithm  is  given  as: 


p(n+t)  =  I(p("))  (25) 

The  algorithm  is  then  given  as  follows: 

Algorithm  1:  Power  Control 

1)  Initialize  p  >  0,  U,  and  V 

2)  Update  p  using  (25)  until  convergence. 

This  algorithm  will  be  used  jointly  with  beamforming 
algorithms  in  future  sections  to  obtain  the  optimal  power 
allocation  for  a  given  set  beamformers. 


B.  Optimal  Beamforming  for  the  Kronecker  Model 
While  for  general  Rayleigh  fading  the  optimal  transmit 
and  receive  beamformers  are  difficult  to  obtain,  the 
problem  in  (12)  can  be  solved  if  the  Kronecker  model  is 
assumed  on  all  the  links.  The  Kronecker  model  assumes 
the  transmit  and  receive  correlations  are  independent 
at  the  link  ends.  This  implies  that  each  transmitter 
has  its  own  transmit  correlation  matrix,  regardless  of 
the  receiver,  and  each  receiver  has  its  own  receive 
correlation  matrix,  regardless  of  the  transmitter.  Thus, 
the  channel  from  transmitter  i  to  receiver  l  is  given  as: 


vec(HK)~CW(0,E;£4®Efl() 

This  assumption  has  some  interesting  implications. 
Note  that  now,  for  the  Kronecker  model,  the  link  gain 
between  receiver  l  and  transmitter  i  is  given  by 


Cu  =  (v*  (8  u ®  s%)(v*  <8>  U i) 

=  (vfs£.vt)(uf : Sr,uj)  =  (vf£TiVi)(uf  Ei^Ui) 

The  last  equality  follows  from  aT  =  a,  where  a  is  a 
scalar  (as  in  the  single-user  case  above).  The  constraint 
equations  in  (12)  then  become 


exp 


(  7 1 _ _ 

V  2  (ViETjVjXufE^uOpt 


X 


n 


<Pi,l=l,...,L 


(26) 


First,  focus  on  the  receive  beamformers.  The  important 
observation  here  is  that  u/  only  appears  in  the  exponen¬ 
tial  noise  term  in  the  product  in  Eqn.  (26),  and  cancels 
out  in  the  other  product  terms,  which  affect  the  interfer¬ 
ence.  Therefore,  the  receive  beamforming  only  helps  to 
mitigate  the  noise,  and  not  the  interference.  With  respect 
to  the  receive  beamforming,  then,  this  case  is  exactly 
analogous  to  the  single  user  MIMO  case,  and  the  optimal 
u;  is  given  by  the  normalized  eigenvector  corresponding 
to  the  largest  eigenvalue  of  E  .  I  =  \ .....  L.  From  this 
solution,  note  that  tEfl,ujopt  =  where  is 
the  maximum  eigenvalue  of  E  . 

For  the  optimal  transmit  beamformers,  define  qi  = 
(vfE Tlvi)pi,l  =  1, ...  jL.  By  noting  that  pi  = 

(v"sT  v/ )  ’  assuming  the  Kronecker  model,  the  following 
equivalent  optimization  problem  to  (12)  is  presented: 


s.t. 


L 

min  ) 

q>0,U,V  ' 

l=i 


wiqi 

(vfETiv,) 


<Pi,l  =  l,...,L 

(27) 


In  the  optimization  problem  in  (27),  the  v/’s  only  appear 
in  the  optimization  function,  and  not  in  the  constraints. 
To  minimize  the  objective  function  in  (27),  then,  the 
denominator  terms  should  be  maximized.  Once  again, 
this  optimization  problem  is  analogous  to  the  single-user 
case,  and  the  optimal  V;’s  are  given  by  the  normalized 
eigenvectors  corresponding  to  the  maximum  eigenvalues 
of  the  Sx;’s.  Since  the  problem  in  (27)  is  equivalent  to 
the  problem  in  (12)  for  the  Kronecker  model,  the  optimal 
transmit  beamformers  will  be  the  same  as  well. 

For  the  special  case  of  i.i.d.  Rayleigh  fading  on  all 
the  links,  e.g.  when  E u  =  I  Vi,*,  it  follows  from  this 
solution  that  any  set  of  transmit  and  receive  beamformers 
is  optimal,  using  the  same  arguments  as  the  single-user 
case.  After  substituting  the  optimal  transmit  and  receive 
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beamformers  into  (12),  Algorithm  1  can  be  run  to  solve 
for  the  optimal  power  allocation,  yielding  the  jointly 
optimal  solution. 

Note  that  due  to  the  limitations  of  the  Kronecker 
model,  the  optimal  transmit  and  receive  beamformers  are 
not  able  to  intelligently  suppress  the  interference  from 
other  users.  This  is  not  the  case  for  the  general  Rayleigh 
fading  case,  as  is  shown  in  the  next  section. 


C.  A  Near-Optimal  Solution  for  General  Rayleigh  Fad¬ 
ing 

The  Kronecker  model  proved  to  be  analytically 
tractable,  but  the  assumptions  made  by  this  model  may 
not  be  adequate  for  many  multiuser  setups-the  trans¬ 
mitter  location  will  affect  both  the  scatterers  that  are 
used  around  both  the  receiver  and  transmitter  in  many 
scenarios.  Both  receive  and  transmit  beamforming  for  the 
general  Rayleigh  fading  case,  however,  are  more  difficult 
to  solve.  As  such,  a  suboptimal  algorithm  is  considered 
in  this  section,  using  bounds  on  the  constraints  to  find  the 
optimal  solution  for  the  ratio  of  the  average  of  the  signal 
to  the  average  of  the  interference  plus  noise.  Looking 
at  this  framework  allows  use  of  modified  versions  of 
previously  developed  techniques  for  algorithms  utilizing 
perfect  channel  knowledge.  This  work  builds  off  of  [7], 
adapting  their  framework  for  the  optimization  problem 
in  (12). 

1 )  Problem  Setup:  First,  consider  the  bounds  given 
in  [21], [22]  for  the  outage  probability,  adapted  to  the 
constraints  in  (12): 

Lemma  4:  Upper  and  lower  bounds  for 
9l( P>  v)  are  given  by 


1+7 1 


X^Z  cliPi  +  2(JNi 


<  gihhP^i,  v) 


<  exp 


CllPl 

7;(E &i  ViPi  + 


CllPl 


The  proof  is  given  in  Appendix  I.  First,  consider 
using  the  expression  for  the  upper  bound  instead  of 
gi^i,  p,uhV)  in  the  constraints  of  (12): 


T 

mm  w  p 

P>0,U,,V; 


s.t. 


\  CliPi  +  \ 

exp  {  7 1 - - -  >  <Pi,l  =  1, 


cupi 


.,L 

(28) 


The  solutions  to  the  optimization  problem  in  (28)  will 
yield  feasible  solutions  to  the  problem  in  (12),  since  the 
upper  bound  will  never  be  less  than  gi(yyi,  p,  u/,  V).  The 
the  constraint  equations  can  then  be  rearranged  to  get 


s.t. 


T 

mm  w  p 

p>0,U,V 
CllPl 


J2i^l  CliPi  +  2a 


12  — 


7 1 


Ni 


log  Pi 


,l  =  l,...,L  (29) 


The  lower  bound  can  also  be  used  for  optimization. 
Using  Lemma  3,  consider  the  optimization  problem  from 
(13).  Substituting  in  the  lower  bound  from  Lemma  4 
for  gi(ji,  p,  u/,  V)  in  the  constraint  equations,  after 
some  manipulation,  results  in  the  following  optimization 
problem: 


T 

max  w  p 

p>0,U,V 

.  cupi  .  li  i  i  T 

S-t.  ^ - ,12  <  -B - =!,■■■, L 

Ei/i  *iPi  +  2<7^  Pi  -  1 

This  problem  can  also  be  stated  as  [7] 


s.t. 


T 

mm  w  p 

p>0,U,V 

CllPl 


Ei^Z  CliPi  + 


I  0-2 
2aN, 


> 


7 1 


Pi-  1 


,1  =  1,...,L  (30) 


The  solutions  to  the  problem  in  (30)  will  not  yield 
strictly  feasible  solutions  to  (12),  but  Algorithm  1  can 
be  applied  to  the  obtained  solution  to  achieve  a  feasible 
solution  for  (12). 

In  most  networks,  the  outage  probability  thresholds 
will  be  small  to  increase  the  likelihood  that  the  links  are 
active.  If  the  outage  probability  threshold  ai  is  small, 
then  pi  =  (1  —  cq)-1  is  larger  than,  but  very  close  to  1. 
Then,  from  the  approximation  log(  1  -\-x)  «  x,  it  follows 
that  log  pi  &  Pi  -  1,  or 


7 1  _  7 1 

log  Pi  ~  Pi  -  1 

Thus,  from  (29)  and  (30),  the  upper  and  lower  bounds 
are  tight  for  low  outage  probabilities. 

Also,  note  that  the  optimization  problems  in  (29)  and 
(30)  take  exactly  the  same  form  as  the  optimization 
problems  from  [7]  (only  certain  constants  are  different). 
Thus  all  the  techniques  from  this  work  can  be  applied 
to  solve  for  the  powers  and  the  transmit  and  receive 
beamformers.  Note  that  while  the  problems  look  similar 
mathematically,  the  quantities  in  the  problem  presented 
in  [7]  and  the  problem  presented  here  vary  significantly- 
the  problem  in  [7]  utilizes  CSI  (thus,  knowledge  of 
{Hjj}),  while  the  problem  presented  here  utilizes  CDI 
(knowledge  of  {S/,}). 


9 


2)  Optimal  Power  Control  and  Beamforming  for  the 
Bounded  Problem:  For  now,  consider  the  optimization 
problem  in  (29).  It  can  be  shown  that  all  the  techniques 
discussed  here  can  also  be  used  for  the  optimization 
problem  in  (30).  For  notational  convenience,  first  define 

1  2 
nl  =  2  aNt 

/  _  ^ 

1  log  fa 

vp;  = _ ML _ 

E i^i  ciiPi  +  m 

Then  (29)  can  be  rewritten  as 

T 

mm  w  p 

p>0,U,V 

s.t.  9i>tf>i,l  =  l,...,L  (31) 


max  d/j 

u; 


CllPl 


cuPi  +  ni 

Pl(yi  ®  U;)HSn(v;*  (8)1 If) 
Ei^ft(v*  ®  U;)  +  ni 


u"#  “'u, 


,1  =  1, . . . ,  L 


(33) 


Note  that  T/  is  only  dependent  on  the  Ith  receive  beam- 
former,  so  the  constraints  can  be  optimized  individually. 
Here,  the  optimal  U  takes  a  slightly  different  form 
than  the  optimal  U  in  [7],  where  the  optimal  receive 
beamformers  are  given  by  the  MVDR  beamformers. 
In  this  problem,  the  above  equation  takes  the  form 
of  the  dominant  eigenvector  for  the  generalized  eigen¬ 
value/eigenvector  problem,  which  is  well-studied  [26]. 
Specifically,  u;  is  the  normalized  vector  that  satisfies  the 
following  equation  for  the  largest  possible  A: 


The  optimization  problem  in  (31)  can  be  rewritten  as  a 
linear  programming  problem  [7],  First,  define  an  L  by  L 
diagonal  matrix  D,  with  elements  =  |M-,  and  a  L  by 
L  matrix  C j  with  cu  as  its  elements  on  the  off-diagonal 
and  0’s  along  the  diagonal.  Using  these  two  matrices, 
also  define  A  =  D  —  C j.  After  some  manipulation,  the 
problem  in  (31)  can  be  rewritten  as: 


$,{n)x  =  A$,(d)x 

Now,  focus  on  the  transmit  beamformers.  The  transmit 
beamformers  can  be  solved  in  an  analogous  fashion  to 
the  receive  beamformers  for  the  dual  optimization  prob¬ 
lem  of  (32).  It  can  be  shown  that  this  dual  optimization 
problem  can  be  formulated  as 


T 

mm  w  p 

P>o,u,v 

s.t.  Ap  >  n  (32) 

The  suboptimal  solution  will  have  a  feasible  solution 
when  U  1 C  /  has  a  spectral  radius  less  than  1  [7],  It  can 
be  shown  that  if  a  feasible  solution  exists,  the  constraint 
inequality  in  (32)  is  satisfied  with  equality,  so  p  is  given 
as: 

p  =  A_1n 

Also,  the  receive  beamformers  can  be  solved  in  this 
setup  for  fixed  transmit  beamformers  and  powers.  To 
minimize  the  objective  function,  'Pi,  l  =  1 .....  T,  should 
be  maximized  for  U.  Doing  this  allows  for  p  (specifi¬ 
cally,  pi  for  the  Ith  constraint)  to  be  minimized  as  much 
as  possible  while  meeting  the  constraints.  Define 

*,(n)  =w(vf  ®I)"Eh(v*®I) 

®  ®  I)  +  nil 

The  following  L  optimization  problems  can  then  be 
written: 


T 

nun  n  q 

q>0,U,V 

s.t.  Arq  >  w  (34) 

Once  again,  using  the  same  arguments  as  before,  the 
constraint  inequality  in  (34)  is  satisfied  with  equality. 
Thus,  the  solution  to  the  optimum  power  allocation  in 
the  dual  domain  is  given  by 

q  =  A  Tw 

Then,  the  optimal  transmit  beamformers  can  be  solved 
for  by  noticing  that  the  Ith  constraint  in  dual  problem 
is  only  dependent  on  the  Ith  transmit  beamformer  (all 
else  held  constant),  for  fixed  receive  beamformers  and  a 
fixed  dual  power  vector  q.  Define  for  the  dual  problem 

*!n)  =«(  I®u;)HS«(I®Ui) 

=  X!  ®  ®  Uj)  +  wil 

i^l 

This  problem  is  analogous  to  the  receive  beamformer 
problem,  so  the  transmit  beamformers  can  be  solved  as 
the  principal  generalized  eigenvector  for  and  . 

From  the  above  arguments,  an  algorithm  for  joint 
power  control  and  transmit/receive  beamforming  can  be 
derived: 
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Algorithm  2:  Joint  Power  Control  and  Beamform¬ 
ing  Algorithm 

1)  Initialize  p  >  0  and  V 

2)  Calculate  the  d?/’s  and  use  them  to  calculate  U 

3)  Calculate  A,  and  then  calculate  q  =  A  7  w 

4)  In  the  dual  domain,  calculate  the  <tvs  and  use  them 
to  calculate  V 

5)  Calculate  A,  and  then  calculate  p  =  A  1  n 

6)  Repeat  steps  2-5  until  convergence. 

This  algorithm  converges  regardless  of  the  initializa¬ 
tion  and  will  generate  feasible  solutions  of  decreasing 
cost  [7],  One  way  to  initialize  the  algorithm  is  to  ignore 
interference  and  use  the  algorithms  for  point-to-point 
communication  given  in  section  III. 

After  Algorithm  2  has  converged.  Algorithm  1  can 
also  be  used  to  achieve  the  optimal  power  control  scheme 
for  the  beamformers  produced  by  this  algorithm.  To 
more  closely  approach  the  optimal  solution,  the  follow¬ 
ing  algorithm  can  be  used: 

Algorithm  3:  Extensive  Joint  Power  Control  and 
Beamforming  Algorithm 

1)  Initialize  p  >  0  and  V 

2)  Run  Algorithm  2  on  (29),  and  run  Algorithm  1 
on  the  result 

3)  Run  Algorithm  2  on  (30),  and  run  Algorithm  1 
on  the  result 

4)  Compare  steps  2  and  3  and  take  the  set  of  beam- 
formers  and  power  vectors  for  initialization  in  step 
1 

5)  Repeat  steps  1-4  until  convergence. 

While  Algorithm  3  will  always  perform  at  least  as  well 
as  Algorithm  2  on  its  own,  the  results  from  the  next 
section  show  this  algorithm  only  provides  marginal  im¬ 
provement.  Algorithm  2  requires  less  run-time,  however, 
since  Algorithm  3  requires  at  least  2  runs  of  Algorithm 
2.  Algorithm  2  is  also  on  the  same  order  of  complexity 
as  the  CSI  algorithm  presented  in  [7],  However,  since 
CDI  changes  much  less  frequently  than  CSI,  algorithms 
based  on  CDI  require  less  feedback  and  need  to  be  run 
much  less  frequently  than  their  CSI  counterparts. 

V.  Results 

A.  Simulation  Parameters 

We  conduct  some  numerical  experiments  to  under¬ 
stand  the  efficacy  of  the  joint  power  control  and  beam¬ 
forming  algorithms  developed  above.  Results  are  given 
for  a  SU-MIMO  setup  and  a  MU-MIMO  setup.  For  the 
covariance  matrices  in  both  cases,  an  angular  spread 
model  with  100  scatterers  and  varying  transmit  and  re¬ 
ceive  angular  spreads  across  the  links  is  considered.  The 
signal  links  are  centered  at  broadside  with  transmit  and 
receive  angular  spreads  varying  from  5  to  20  degrees. 


and  the  interfering  links  are  centered  at  incident  angles 
with  transmit  and  receive  angular  spreads  varying  from 
10  to  40  degrees.  The  parameters  were  selected  with 
respect  to  typical  scenarios  in  the  3GPP  model  [27], 
All  transmitters  and  receivers  have  4  antennas,  and  the 
signal-to-noise  ratio  (SNR)  is  held  constant  at  10  dB. 
Equal  SINR  thresholds  and  equal  outage  constraints  is 
considered  for  all  users  in  the  system. 

For  the  MU-MIMO  setup,  a  network  with  3  trans¬ 
mit/receive  links  and  a  weighting  vector  w  =  [10, 1, 1]T 
is  considered. 

B.  Single-User  MIMO 

In  Experiment  1,  the  single-user  MIMO  case  is  con¬ 
sidered.  Here,  an  algorithm  using  perfect  channel  knowl¬ 
edge  is  compared  to  the  algorithm  presented  in  Section 
III,  which  uses  the  covariance  information  only.  The 
algorithm  using  perfect  channel  knowledge  uses  the 
principal  left  and  right  singular  vectors  for  receiver  and 
transmitter  beamforming,  respectively,  and  calculates 
the  minimum  transmit  power  required  to  achieve  the 
threshold  (this  is  an  optimal  scheme).  The  outage  value 
at  which  the  average  power  required  using  CSI  matches 
the  power  required  using  CDI  is  a  function  of  the  inverse 
of  the  expected  dominant  singular  value  of  H,  and  can 
be  empirically  calculated  for  a  given  covariance  matrix. 
See  Fig.  1. 

The  plot  yields  some  interesting  results:  the  average 
transmit  power  required  for  CDI  is  comparable  to  the 
average  power  required  for  CSI  for  the  same  threshold, 
with  the  value  being  lower  at  20%  outage.  To  calculate 
the  CDI  powers,  only  the  covariance  information  is 
required,  so  it  needs  to  be  updated  only  when  the 
channel  statistics  change.  On  average,  2.5dB  less  power 
is  required  to  transmit  using  CDI  at  20%  outage  as 
compared  to  CSI.  The  tradeoff  is  that  the  covariance 
algorithm  will  allow  the  link  to  fail  with  the  specified 
outage  probability.  However,  this  may  be  a  desirable 
tradeoff  for  many  network  applications  since  it  prevents 
the  transmitter  from  attempting  to  increase  the  power 
to  compensate  for  poor  channel  conditions,  thereby  de¬ 
creasing  the  interference  to  other  nodes  in  the  network. 

C.  Mu Iti-User  MIMO 

For  the  multiuser  case,  the  algorithm  from  [7]  is  used 
to  calculate  the  optimal  transmit/receive  beamformers 
and  powers  at  each  channel  instance  (using  perfect 
channel  knowledge).  In  Experiment  2,  the  average  value 
of  the  objective  function  using  the  weighting  vector 
above  is  compared  against  the  suboptimal  Algorithm  2 
for  varying  outage  constraints.  See  Fig.  2. 

This  experiment  shows  similar  results  to  the  single- 
user  case,  except  the  solutions  generated  by  Algorithm 
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2  become  infeasible  (when  the  curves  go  to  oo).  The 
SINR  thresholds  here  are  reasonable,  as  many  cellular 
networks  today  are  able  to  operate  in  the  low  SINR 
regime  [27],  The  lower  the  outage,  the  lower  the  SINR 
threshold  has  to  be  for  the  outage  to  become  infeasible. 
Once  again,  the  value  of  the  average  objective  function 
is  comparable  to  using  true  channel  knowledge  for  low 
enough  SINR  thresholds.  This  means  that  knowing  only 
covariance  information,  if  the  system  is  allowed  to  fall 
into  outage  with  some  probability,  transmission  can  be 
performed  on  average  with  comparable  power  than  with 
knowing  perfect  channel  information  all  the  time.  Using 
Algorithm  2,  there  is  about  a  ldB  savings  in  the  average 
cost  function  at  an  SINR  threshold  of  -4dB  at  20% 
outage  compared  to  the  CSI  scheme. 

For  more  insight  into  the  power  consumption  using 
perfect  channel  knowledge  and  covariance  information, 
also  observe  the  CDF  of  the  weighted  sum  power  for 
true  channel  knowledge  at  an  SINR  threshold  of  0 
dB  in  Fig.  3.  This  plot  shows  what  percent  of  the 
time  the  objective  function  is  larger  for  perfect  channel 
knowledge  as  compared  to  the  outage  curves.  The  outage 
curves  are  just  step  functions,  since  for  fixed  covariance 
information,  they  will  always  have  the  same  value  for 
their  objective  function.  Since  the  CDF  curves  for  perfect 
channel  knowledge  and  20%  outage  intersect  at  roughly 
80%,  the  data  shows  that  roughly  20%  of  the  time, 
using  perfect  channel  knowledge,  the  objective  function 
is  higher  than  Algorithm  2  for  20%  outage.  Also  observe 
that  the  median  (where  the  CDF  takes  a  value  of  0.5) 
is  approximately  -23  dB  on  the  CDF  plot.  The  mean 
from  Fig.  2,  however,  is  at  approximately  -17  dB,  a 
difference  of  6  dB,  or  a  factor  of  about  4.  The  large 
disparity  between  the  median  and  the  mean  suggest 
that  the  data  for  perfect  channel  knowledge  is  skewed 
when  the  channel  is  bad.  Thus,  using  Algorithm  2 
significantly  reduces  the  peak-to-average  power  ratio  for 
the  transmitters  (the  peak-to-average  power  ratio  is  1 
using  Algorithm  2  since  the  power  allocation  does  not 
change  with  the  channel). 

Experiment  3  shows  the  the  difference  between  using 
Algorithm  2  and  Algorithm  3.  See  Fig.  4.  The  curves 
show  that  the  objective  function  is  nearly  identical  with 
and  without  optimal  power  control,  and  Algorithm  3 
only  begins  to  give  any  discernible  gains  when  they 
approach  infeasible  SINR  thresholds.  This  is  consistent 
with  the  analysis  of  the  tightness  of  the  bounds  in  Section 
IV-C. 

Experiment  4  tests  the  convergence  speed  of  Algo¬ 
rithm  2.  See  Fig.  5.  For  the  setup  in  the  experiment,  both 
the  outage  values  as  well  as  the  objective  function  values 
converge  within  5  iterations  to  their  fixed  point  solutions. 
This  confirms  that  this  algorithm  is  nearly  identical 


in  computational  complexity  to  the  algorithm  used  for 
perfect  channel  knowledge  [7],  The  big  difference  here, 
however,  is  that  Algorithm  2  needs  to  be  run  much 
less  frequently  than  the  algorithm  using  perfect  channel 
knowledge. 


VI.  Conclusion 

In  this  work,  a  framework  has  been  presented  for 
analyzing  MIMO  beamforming  networks  when  only  the 
channel  statistics  are  known,  and  algorithms  have  been 
derived  that  give  good  performance.  Optimal  algorithms 
have  also  been  presented  for  the  single  user  MIMO  case, 
as  well  as  for  the  multiuser  MIMO  Kronecker  case. 
A  suboptimal  algorithm  using  duality  for  joint  power 
control  and  beamforming  for  the  general  Rayleigh  fading 
case  was  also  derived.  While  finding  the  optimal  beam- 
formers  for  this  problem  is  an  open  area  of  research, 
given  the  tightness  of  the  bounds  and  the  negligible 
improvement  of  using  the  optimal  power  control  algo¬ 
rithm  on  top  of  the  presented  suboptimal  algorithm,  the 
gains  in  finding  the  optimal  beamformers  will  likely  be 
negligible  for  reasonable  outage  constraints. 

When  comparing  the  algorithms  presented  here  to 
algorithms  where  the  channel  is  perfectly  known,  it 
has  been  shown  that  using  covariance  information  can 
give  gains  in  the  cost  function  of  interest.  Thus,  the 
solution  presented  here  has  two  substantial  advantages: 
the  computational  complexity  for  the  algorithm  and  feed¬ 
back  information  required  for  covariance  information  is 
drastically  reduced  since  both  must  be  done  only  once 
for  valid  covariance  information,  and  transmission  power 
can  be  saved  in  the  network  if  a  sufficient  level  of 
outage  is  acceptable  on  the  links.  Also  note  here  that 
the  algorithms  presented  here  can  be  extended  to  decen¬ 
tralized  algorithms  as  proposed  in  [7],  among  others,  in  a 
straightforward  manner,  though  this  is  still  an  open  area 
of  research.  Knowledge  of  covariance  information  has 
many  benefits,  and  it  can  be  an  effective  way  to  reduce 
feedback  in  MIMO  networks. 


Appendix  I 
Proofs  of  Lemmas 

Proof  of  Lemma  1:  First,  consider  Z  =  uffHv. 
Z  is  a  linear  combination  of  0-mean  complex  circular 
Gaussians,  so  Z  is  0-mean  complex  circular  Gaussian. 
The  next  step  is  to  find  the  variance.  To  do  this,  perform 
the  following  calculations: 
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var(Z)  =  E[Z  •  Z*]  =  £[(uffHv)(uffHv)*]  (35) 

=  £[vec(uHHv)vec(uffHv)H]  (36) 

=  E[(vT  ®  uif)vec(H)vec(H)ff (v*  ®  u)] 

(37) 

=  (v*  <g>u)ff£(v*  <g>u)  (38) 

The  key  observation  here  is  that  equality  (36)  follows 
from  observing  that  x  =  vec(x)  if  a:  is  a  scalar. 
Furthermore,  for  scalars,  the  Hermitian  and  conjugate 
operations  are  interchangeable.  Equality  (37)  follows 
from  the  property  that  vec(AXB)  =  (BT  (g>  A)vec(X). 
Finally,  equality  (38)  follows  from  separating  out  the 
deterministic  factors  from  the  expectation  operator  and 
noting  that  by  definition,  .E[vec(H)vec(H)H]  =  S.  This 
gives  an  expression  for  var(Z). 

The  norm  squared  of  a  unit  variance  complex  circular 
Gaussian  has  a  or  exponential,  distribution.  Then 
\Z\2,  normalized  by  vcir(Z ),  gives 

\Z\2  _  |ugHv|2  2 

var(Z)  (v*  0  u)HS(v*  ®  u)  ~  *2 

□ 

Proof  of  Lemma  2:  First,  use  contradiction  to  prove 
2).  Note  that  p,  u;,  V)  in  (12)  is  a  decreasing 
function  of  pi.  Now,  assume  that  for  the  Ith  constraint, 
gii.'lh  Pi  uh  V)  is  strictly  less  than  f3i  in  the  optimal 
solution,  and  all  the  other  inequalities  are  satisfied.  This 
means  pi  can  be  reduced  appropriately  to  get  equality 
and  reduce  the  value  of  the  objective  function.  Further¬ 
more,  since  g-f  y-i,  p,  ui;  V)  is  an  increasing  function  in 
Pi,i  l,  when  pi  is  reduced  all  of  the  other  constraints 
are  still  satisfied.  Thus,  having  the  Ith  constraint  not 
satisfy  equality  contradicts  optimality,  so  all  constraints 
are  satisfied  with  equality. 

Also  note  here  that  if  only  the  Ith  constraint  holds  with 
strict  inequality  and  the  rest  of  the  constraints  hold  with 
equality,  when  pi  is  reduced,  all  of  the  other  constraints 
become  inequalities  (since  gt  (7* ,  p,  U;,  V)  is  an  function 
in  pi ,  i  f  l ).  Thus,  the  rest  of  the  pf  s,  i  ^  l,  can  be 
reduced,  which  will  lead  to  a  further  reduction  of  pi, 
until  a  stable  solution  is  reached  that  gives  equality  on 
all  the  constraints. 

Contradiction  can  also  be  used  to  prove  1).  Assume 
there  exists  two  distinct  optimal  solutions  that  minimize 
the  cost  function,  p  and  p*.  Then,  define  k  and  a  such 
that  a  =  ma ttipi/p*  and  k  =  argma xpi/pf  Since  p 
and  p*  are  distinct  and  both  minimize  the  cost  function 
(so  wTp  =  w7  p*),  at  least  one  element  in  p  is  greater 
than  its  corresponding  element  in  p*  (it  also  follows  that 
at  least  one  element  in  p  is  less  than  its  corresponding 
element  in  p*).  Therefore,  a  >  1.  Then,  using  2)  from 


this  lemma,  consider  the  kth  constraint,  satisfied  with 
equality: 


Ik  Nk. 


fa  =  gkhk,p*,uk,  v)  =  e 2  okkpi 


i^k 


Ik 

>  e  2  CfcfcPfca 


i^k 


2 b.  Nk 

=  e  2  °kkPl 


ftk  V  ckkPka) 


2 k  Nk 
=  e  2  ckkPk 


i^k 


n  i+ 


■7k~ 


CkiPj  a 

CkkPk 


ik  Nk 
>  e  2  ckkPk 


-LX  V  CkkPk  ) 


i^k  x 

=  9k(lk,p,  Ufc,  V)  =  /3fc 


CkiPt  \ 
CkkPk  ) 


(39) 


(40) 


(41) 

(42) 

(43) 

(44) 


The  inequality  in  (42)  comes  from  the  fact  that  p*a  > 
Pi  for  all  i  because  of  how  a  was  selected,  and  that 
9i('liiPiuh'V)  is  an  increasing  function  in  pi,i  l. 
Because  of  the  inequality  in  (40),  the  above  steps  show 
a  contradiction,  so  p*  must  be  unique. 

There  exists  a  unique  solution  for  L  equations,  so 
the  constraints  completely  determine  p*,  regardless  of 
w.  Thus,  3)  follows  immediately  from  1)  and  2),  which 
concludes  the  proof. 


□ 


Proof  of  Lemma  3:  Similar  arguments  to  the  ones  used 
to  prove  Lemma  2  can  be  used  for  this  lemma  as  well. 
All  the  constraint  equations  will  hold  with  equality  by 
the  same  argument  as  given  previously,  and  the  solution 
is  unique,  being  solely  determined  by  the  constraints. 

□ 


Proof  of  Lemma  4:  First,  consider  the  lower  bound. 
A  well-known  bound  is  log(l  +  x)  <  x,  or  1  +  x  <  ex . 
Applying  this  bound  results  in 


exp 


ZcuPi 


>  1  + 


2  cupi 


Then  consider 


9l{lhPi  V) 


CliPi  \ 
CllPl  ) 


(45) 
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The  expression  resulting  in  (45)  can  be  expanded: 


1  + 


=  1 


>  1 


7*4,  | 

2 cupi  J 

Hid 

i^l 

/ 

CliPi 

+  7/ 

CllPl 

+  7t4A 

2 cupi  J  \ 

1  +  ^7/ 

i^l 

7/4, 

,  V 

^  CliPi  ! 

2  cupi 

i^l 

CllPl 

7/4, 

,  V 

CliPi 

2  cupi 

i^l 

ll 

cupi 

CliPi 

CllPl 


(46) 


\ '  CuPi 

(47) 


i^l 


CllPl 


=  1+  7 i 


cliPi  “I"  2<JNi 


(48) 

(49) 


CllPl 

In  Eqn.  (46),  the  inequality  comes  from  expanding  the 
product  on  the  left  hand  side  of  the  equation  and  only 
keeping  the  two  terms  shown  (the  terms  excluded  after 
expanding  are  positive,  resulting  in  the  inequality).  The 
inequality  in  Eqn.  (48)  once  again  comes  from  excluding 
a  positive  term  in  the  sum.  This  gives  the  lower  bound. 

For  the  upper  bound,  once  again  consider  the  bound 
1  +  x  <  ex.  Applying  this  to  the  product  terms  gives 


2L  -I — r  I 

#(7i,p,uz,  V)  =  e  2  °u*i  1  +  7* 


i^l 


CliPi 

\ - 

'  cupi 


21  Ni 


<  e  2  c“pi  e7i  cnpi 


i^l 


=  exp  \  7 1  - 

This  concludes  the  proof. 


CliPi  +  2° 


1^2 


2UJV, 


CllPl 


□ 


Appendix  II 

Derivation  of  Standard  Interference 
Function  Properties 

Standard  interference  functions  were  first  defined  in 
[28],  In  this  work,  the  authors  showed  that  this  class  of 
functions  can  be  used  to  do  optimal  power  control,  and  it 
has  nice  convergence  properties.  The  definition  is  given 
as  follows: 

Definition  1:  An  function  F(p)  is  a  standard  interfer¬ 
ence  function  if  the  following  properties  are  satisfied: 

1)  Positivity:  F(p)  >  0 

2)  Monotonicity:  If  p  >  p',  then  F(p)  >  F(p') 

3)  Scalability:  For  all  k  >  1,  «F(p)  >  F(/cp) 

Theorem  2:  Define 


h{v) 


9+4, 

2  cu  log  fa 


Pi 

log  Pi 


Elog 

i^l 


+  7 1 


CliPi  \ 
cupi ) 


(50) 


Then,  the  function  I(p)  =  [/i(p), . . . ,  /l(p)]  is  a 
standard  interference  function. 


Proof:  A  similar  function  and  proof  is  given  in  [22], 
Start  with  the  first  property,  positivity.  This  can  be  proved 
by  showing  all  the  terms  in  Eqn.  (50)  are  nonnegative, 
and  at  least  one  of  them  is  positive.  7 1  and  <j2Nj  are  >  0 
by  how  they  are  defined.  From  the  constraints,  pi  >  0 
(the  power  can  never  be  negative).  Also,  the  logarithm  of 
any  number  greater  than  1  is  positive,  pi  =  (1  —  a;)-1, 
and  0  <  ai  <  1,  so  /?,  >  1.  Lastly,  cu  is  positive  since 
cu  =  (v* 0U;)  >  0  since  is  a  positive 
semidefinite  matrix.  It  is  also  safe  to  assume  that  u,  and 
v,  will  be  selected  such  that  cu  >  0,  since  if  cu  =  0 
the  expression  for  outage  probability  becomes  undefined. 
Therefore,  all  the  values  in  the  first  term  of  the  sum  in 
(50)  are  strictly  greater  than  0  and  the  second  term  in 
the  sum  is  >  0,  so  therefore  I(p)  >  0. 

Monotonicity  can  be  proved  by  showing 
E^;log(l  +  7;fift)  is  monotonic  in  p 
(since  the  other  term  in  Eqn.  (50)  is  a  constant  with 
respect  to  p).  First,  focus  on  only  one  term  in  this  sum, 
since  a  sum  of  monotonically  increasing  functions  will 
also  be  monotonically  increasing.  Then,  consider  the 
function 


h{pi,Pi) 


Pi 

log  Pi 


log 


(51) 


If  pi  is  held  constant,  then  h(pi,Pi)  is  monotonically 
increasing  in  pi,  since  clog(l  +  kx)  is  monotonically 
increasing  for  a  fixed  c  >  0  and  k  >  0.  If  p,  is  now  held 
constant,  h(pi,pi )  takes  the  form 


2;  log 


x(log(x  +  k)  -  log(x)) 


(52) 


The  function  h  (pi  -  p,  )  is  also  monotonically  increasing 
in  pi  for  k  >  0  since  Eqn.  (52)  always  has  a  positive 
derivative  for  x  >  0.  Therefore,  I(p)  is  monotonically 
increasing. 

To  prove  scalability,  consider  a  constant  k  >  1.  Then 
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h(n P) = 


WNl 


WNl 


<  F- 


7  l°Nt 


KPl 


2 cu  log  (3 1  log  pi 


lQg  (  1 


lv 


Pi 


2 cu  log  Pi  log  pi 


Y  lo§  (  1 


Pi 


2  cu  log  Pi 

=  kIi(  p) 

This  concludes  the  proof. 
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Single-User  4x4  MIMO 


Fig.  1.  Single  User  Case,  SINR  threshold  vs  Power 


Multiuser  4x4  MIMO,  L=3:  CDI  vs  CSI  Algorithms 


Fig.  2.  Multiuser  Case,  SINR  threshold  vs  Weighted  Sum  Power 


Multiuser  MIMO,  L=3,  SINR  Threshold  =  OdB: 


Cumulative  Distribution  Function  of  Weighted  Sum  Power 


Fig.  3.  Multiuser  Case,  CDF  of  Weighted  Sum  Power 
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Multi-User  MIMO,  L=3:  Algorithm  2  vs  Algorithm  3 


Fig.  4.  Multiuser  Case,  Algorithm  2  vs  Algorithm  3 


Multi-User  4x4  MIMO,  L=3,  SNR  Threshold=-4dB: 
Weighted  Sum  Power  vs.  Iterations 


Multi-User  4x4  MIMO,  L=3,  SNR  Threshold=-4dB: 
Outage  Probability  vs.  Iterations 
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Fig.  5.  Multiuser  Case,  Algorithm  2  Convergence 


